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1 Introduction 

We present here an effectively nonlinear, non-Lagrangian system of PDE which follow from an invariant 
differential-form condition 

d£ = <S>dX£ (1) 

for a 2-spinor and a Mat(2, C)-valued gauge $>(x) fields. Condition (|l]) has been proposed in 
0] in the general framework of algebrodynamical approach to field theory []. In (Q) d is an operator 
of external differentiation, and by X a hermitian 2x2 matrix of the Minkowski space-time coordinates 
is denoted. System (|T|) originates as a particular case of the differentiability conditions for functions 
of biquaternionic (B-) variable (we recall that the algebra of complex quaternions IB is isomorphic to 
the full algebra of 2 x 2 complex matrices, IB = Mat{2, C)). Geometrically, it defines a spinor field 
covariantly constant with respect to an effective IB-valued connection 1-form = &dX 0. 

Remarkably, system (|]) is overdetermined, so that both the 2-spinor and the gauge fields are to 
be found from it in a self-consistent manner. The integrability conditions for (||) take the form 

dd£ = = (d$ - $dX&) A dX£ = R£ (2) 

1 ln algebrodynamics one claims to derive field equations from only the intrinsic properties of fundamental mathemat- 
ical structures, in particular from exeptional quaternion-type algebras. 

2 In a 4-vector representation this connection 1-form gives rize to an effective Weyl-Cartan geometry with the Weyl 
nonmetricity vector and the pseudotrace of the torsion tensor being proportional to each other, see section 5. 



and, as shown below, lead to self-duality of the curvature 2-form R of the connection f2. Consequently, 
the source- free Maxwell and Yang- Mills equations are satisfied identically on the solutions to (|l|), for 
the trace and the trace-free parts of the curvature R respectively. For this reason, in what follows 
system ([!]) is referred to as the generating system of equations (GSE). 

Moreover, the 4-eikonal and the linear wave equation do hold both for each spinor component 
{£,a}, A = 0,1 (the latter only for the quotient of components). Also, a null geodesic congruence 
which may be constructed via its tangent 4-vector k^ = £,a£,A' turns out to be shear-free on the 
solutions of (|T|) and, therefore, the Kerr theorem [14] may be applied to the GSE. On this track, 



general analytical solution to ([[]) may be obtained in terms of twistor variables T = {^,iX^} and in 
an implicit algebraic form 

n( c )(£,iX£) = 0, (3) 

where {lK c )}, C = 0, 1 are two arbitrary and independent holomorphic functions of 4 twistor com- 
ponents. If some n(0 are chosen, a solution £a(x) can be extracted from the system of two algebraic 
equations (H|), and all the analytical solutions to the GSE can be obtained in this way. 

Expression ([|) generalizes the Kerr theorem which deals with only one homogeneous function 
II of 3 projective twistor components. From the full twistor structure of general solution (||) gauge 
symmetry of the GSE (|lj) of a rather peculiar type does follow, the so called "weak" one, when the 
gauge parameter is allowed to depend on coordinates only implicitly, i.e. via the spinor £( ) being 
transformed and its twistor counterpart ^£(0)- 

Now, in a usual way a Riemann metric g of the Kerr-Shild type may be defined through the 1-form 
of the congruence k = k^dx^ as 

g = f] + hk(g)k, (4) 

with r/ being the flat Minkowski metric. Well known is the fact that for a shear-free congruence 
the metric (||) often satisfies the Einstein or Einstein-Maxwell (electro)vacuum equations under an 
appropriate choice of the "potential" function h(x) |]33|| . 

Thus, at least the source-free Maxwell, Yang-Mills and the metric (effective gravitational) fields 
can be naturally assosiated with the solutions of the GSE ([l|). Wonderfully, singularities of strengths 
(of curvatures) of all these fields coincide in space and time, being determined by a single condition 

det 1-3^-11 = 0, (5) 

which specifies the points where the solution to (||) becomes not unique. Geometrically, these points 
constitute a caustic of a correspondent shear-free light bundle. We'll see that the singular set, as 
defined by (||), may be point-like, string-like and even a two-dimensional surface. In case when this 
set is bounded in 3-space we can identify it with a particle-like object whose time evolution is fully 
governed by the field equations (|l|) and may be obtained from the algebraic system @,(|5|). 

These singular objects manifest, at a purely classical level, some properties of the real quantum 
particles; in particular, the value of electric charge q is either zero or a whole multiple of the charge 
of the fundamental static solution; the latter corresponds just to the Kerr shear-free congruence with 
the ring-like singularity. The assosiated metric and electromagnetic field are identical to those of the 
Kerr-Newman solution in GTR, exept for the restrictions on the admissible value of electric charge 
qo = ±1/4 so that the quantity q$ can be naturally identified with the elementary charge. In 
account of the well known property of the Kerr-Newman ansatz to have anomalous Dirac value of the 

3 In the dimensionless units we use; the numerical value itself is of no significance. 
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gyromagnetic ratio we conclude that the fundamental static solution to the GSE reproduces all of 
the quantum numbers of the real electron. 

Throughout the whole paper, we make no assumptions of either physical or mathematical char- 
acter. Nor we introduce some additional dynamical equations for the fields involved, or use some 
Lagrangian structure. What we actually do is we study the intrinsic mathematical properties of the 
"generalized Cauchy-Riemann equations" as they are (precisely, a particular subclass of their solu- 
tions, represented by GSE ([[])). Attempts to use such structures in theoretical physics have been 
repeatedly undertaken, especially in the works of F.Giirsey and his co-workers |3], |j where the self- 
dual gauge and the chiral structures have been reformulated on the basis of quaternionic calculus. 
Most of considerations, however, made use of the only recognized version of quaternionic analiticity 
constructed by R.Fueter 13]. Unfortunately, his approach leads to trivial linear C-like structure of 
the differentiability conditions ignoring the noncommutativity of the quaternion-type algebras and 
considerably limiting possible physical applications of the theory. 

In the paper we don't intend to discuss the problem of noncommutative analysis in detail, for 
this referring the reader e.g. to the profound paper of A.Sudbery |J. In the appendix we only 
frmulate the conditions of B- differentiability in the version earlier proposed in [|], ^, |2|], ^] and 
derive the GSE ([!]) as a particular case of them. To our knowledge, it is the only approach in which 
just the nonlinear structure of the generalized CR-equations naturally emerges, as a consequence 
of noncommutativity of the basic algebra (the situation well known for the gauge groups but quite 
unfamiliar in the analysis). In view of this, one should not be surprised that the field dynamics 
induced by the nonlinear differentiability conditions appears to be quite nontrivial and, perhaps, has 
something to do with real dynamics of particles. Anyway, the main goal of this article is to develop an 
unexpectedly rich and self-consistent "virtual" physics which is contained in a concise and invariant 
structure of the B-differentiability conditions (||) themselves, despite any additional phenomenological 
assumptions. 

The structures which arise in this context have some features in common with the I.Robinson's 
theorem Q relating shear-free congruences to (null) electromagnetic fields; with the approach of 
R.Penrose to the theory of spin-3/2 field relating twistor space to the charge space of electromagnetic 
field asssociated jjj and, especially, with spinor connections introduced by K.P.Tod Q. From the 
generic viewpoint, our approach is adjacent to that of Yu.Manin and G.Henkin [|| in which the 
nonlinear Dirac-Maxwell system is reduced to the linear Cauchy-Riemann equations in twistor space. 
On the whole, however, our approach does not appeal to any others and new physical concepts arising 
therein are not made up or artificial but quite inevitable in a consecutive mathematical scheme used. 

Let us sketch out the organization of the article. In the second section we establish the relativistic 
invariance of the GSE and obtain the 4-eikonal equation for the components of the 2-spinor. In section 
3 we demonstrate the functional dependence of the components of a twistor assosiated with the GSE 
and thus reduce the latter to the system of algebraic equations (^). The fourth section is devoted to the 
analysis of the gauge symmetry of the GSE and to its specific relation to the projective transformations 
in the associated twistor space. In section 5 we study the self-dual structure of the GSE which follows 
from its integrability conditions and guarantees that the source-free Maxwell and Yang-Mills type 
equations are satisfied on the solutions of the GSE. In section 6 we establish the relations of GSE to 
the shear-free geodesic congruences and to the Riemann metrics of the Kerr-Schild type. The review 
of the " particle- like" solutions of GSE and the property of electric charge quantization are presented 
in section 7. In the last section 8 we conclude by the discussion of general problems and perspectives 
of the theory. In the appendix we formulate the conditions of differentiability for A-valued functions 
of a noncommutative algebraic variable Z £ A and, for A = B, derive the GSE(l) as a particular case 
of these conditions. 

Throughout the paper (exept in the appendix) the standard two-component spinor formalism is 
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used. Upper case Latin indices range and sum over zero and one, and are raised and lowered by the 
symplectic spinors e AB and cab respectively, as in 14]. By V aa' the usual spinor derivative operator 
in Minkowski space is denoted. 



2 Spinor structure and the eikonal equation 

In the matrix representation of the biquaternion algebra B used in (|l|) we regard £ as a column (2x1 
matrix over C), $ - as a 2 x 2 complex matrix of general type, and X - as a hermitian matrix 
representing the coordinates of Minkowski space-time. System (|l|) is evidently invariant under the 
global Lorentz transformations of coordinates 

X -> LXL + , £ -> Z + £, $ -» L + $L, § -> L$L+, (6) 

where L £ SL(2, C) and <£> denotes the matrix adjoint to <£: = (det (we admit here the 
possibility for det<I> to be zero). We see that the field £(ic) transforms as a (conjugate) 2-spinor 
whereas the components of the field $>(x) constitute a complex Lorentz 4- vector (we'll identify it later 
with the electromagnetic 4-potential). 

In view of the form of transformations (|6|), functions £, <3? are mapped respectively to matrices £a' 
and $ b'Ai whereas dX - to a hermitian matrix dX AA , with (un)primed indices A, ...A' , ... = 0, 1; 0', 1' 
having usual 2-spinor sense. In matrix notation, (||) has then the form 

<%B> = <Z>B'AdX AA 'U', (7) 

equivalent to a system of eight PDE 

Va4'£b' = $b'a€a'- (8) 

Throughout the paper, we assume for both spinor components £a' to be nonzero in the region of 
spacetime considered (otherwise, the solution to (|l|) can be proved to be degenerate in some sense, with 
zero electromagnetic and other fields associated). We also consider all the functions {£a'( x )i ^b'a(x)} 
to be analytical everywhere exept maybe at some subset of zero measure where they are allowed to 
turn to infinity. 

Some properties of the solution £(x) can be inferred directly from (|8|). Using the orthogonality 
identity £ A 'U' = one easily find 

V aa '&V A a<b' = 0, (9) 
which in particular implies the eikonal equation for any function \(£a') °f spinor components 

v aa '\(U'Waa'HU') = o. (io) 

Returning to system (pi), multiplying it by £ A and again taking into account £ a '£,a' = we reduce 
it to the following form fl27] : 

t A 'V A A>tB> = 0, (11) 

where the 4- vector field $b'A has been dispensed with. The latter may be recovered by using (g) with 
equal indices A' = B' 

&A'A = Vaa'Ih&i'; ( 12 ) 

4 We recall that both components of the spinor considered are assumed to be nonzero: otherwise all of the strengths 
(curvatures) vanish identically. 
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certainly, no summation over A' is assumed in (p^). 

Conversely, from (|Tl| ) in view of a well-known property of 2-spinors it follows 



(13) 



with some spintensor (fB'A- I n compare of (13) and (|8|) we conclude with the equivalence of (11) to 
the original spinor system (||) and, therefore, - to the GSE (|]). 



3 Twistor structure and general solution of GSE 

We turn now to the solutions of ( JTT| ) or of GSE (^) equivalently. Remarkably, they happen to be 
completely determined by a twistor structure which can be naturally assosiated with the system 
considered. To demonstrate this, let us introduce another 2-spinor t a related to £a' via the Klein- 
Penrose correspondence 

r A = X AA 'U>. (14) 

Then a pair of 2-spinors T a = (£,a'i tA )i a = 0,l,(y,l', constitute a (null) twistor incident with a 
(real) Minkowski space-time point represented by X AA . 

According to definition ( |l4|) the wedge product of the differentials d^A' an d dr A may be formed as 

dr A A d£ B > = X AA 'd£ A > A d£ B , + U'dX AA ' A d£ B , . (15) 



From (15) a rather obvious property about (nontrivial) twistors immediately follows: at least some two 
components of a generic twistor T a should be functionally independent (as functions of coordinates 
X AA '). Indeed, assuming conversely for all exterior products of the differentials d^A' an d dr A to vanish 
due to their functional dependence, we get from (|l5|): either both £a' = or VaA'Cb' = 0, but in the 
second case both £a' are constant and the two components of the spinor t a ([14|) are then evidently 
independent. 

If we subject now the spinor £4/ to the dynamical system (O) the remaining two components of 



the twistor should depend on the first independent two. Precisely, we intend to prove the following: 



Preposition. 1 Iff ^A' is a solution of jjj[ ) then the corresponding twistor T a has two and only two 
functionally independent components. 

Going to differentials in ( |i"I|) and using (0), we come to equations for dr A , which together with (|7|) 
themselves constitute a system of four equations 

(%B> = $B>AW A , (16) 

dr B =X BB '^> B/A w A + w B , (17) 

where the 1-forms w A = dX AA £a' have been introduced. Since the differentials of twistor components 
are linear functions of the two 1-forms w A only, it becomes obvious that the exterior product of any 
three is zero, resulting in the desired functional dependence. 

The same conclusion could be reached in a slightly different way. Since two twistor components 
are certainly functionally independent, we can always find two equations in (|l6l) , ([l7|) which allow 
to express w A through the differentials of the corresponding twistor components. Substituting the 
resulting expressions in remaining two equations, we end up with two relations each containing three 

5 For our purposes we may ignore imaginary unit i in definition (|l^). 
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of the four differentials d^A' and dr A . These relations imply the sought-for functional dependence 
between any three twistor components, a fact that can be expressed in a more symmetric form [ p7f ] 

n (C)(T a ) = Ii {c \i A ,,T B ) = 0, (18) 

where {n^}, C = 0, 1 are two arbitrary but independent holomorphic functions of four complex 
variables. 

Conversely, the algebraic system ( |l8|) implicitly determines ^a' and it is easy to check, by differ- 
entiating (18) and multiplying by £ A , that satisfies the system 



B' 



r vaa'^ == o, (19) 



which is equivalent to (11) except at some singular points (see below). Successively resolving system 
(18) at each space-time point X AA with respect to £4/ and substituting the resulting solution in 
(12) to find the corresponding "potentials" &b'A we obtain a solution to the GSE starting from the 
algebraic constrains fll8|). This furnishes the proof of preposition. 1. ■ 

Thus, algebraic system of equations (|l8| ) implicitly determines the general (analytical) solution 
{£,A'i&b'a} of the GSE. Points where the equations ( |l8| ) have multiple roots, i.e. cannot be in a 
unique way resolved for £a' satisfy according to (|l9| ) the equation 

,.dn (c) ,, , N 

det ||— - — 1|=0. (20) 
d^A* 

These points constitute a singular set for electromagnetic field which in the next section will be 



assosiated with the quantities &a'A- Together with (IS) the last equation allows us to determine the 
shape and the time evolution of singularities (see below, section 7). 

Geometrically, the algebraic system ( |i~8| ) defines a two-dimensional complex surface in the twistor 
space C 4 (precisely, in the subspace of null twistors). Points of intersection of this surface with two- 
dimensional planes formed by (null) twistors (|l4j) represent the solution £a' to the GSE (multivalued 



in general) for each fixed space-time point X . Singularities are then the pre-images (in M) of the 



points of twistor space at which the planes (14) are tangent to the surfaces fllq), so that the singular 
set will be the same for all of the branches of a multivalued solution. Note that we ignore here the 
generally considered projective structure of the twistors which, in the framework of this approach, has 
some peculiarities and is related to an exotic gauge symmetry of GSE. These issues are discussed in 
the next section. 



4 Projective transformations of twistors and "weak" gauge symme- 
try of the GSE 

For an appropriate electrodynamical interpretation we need to establish the gauge invariance of the 
GSE, which will be dealt with in this section. Specifically, we shall study the symmetry of (JsJ) under 
transformations 

U>^?a>=<*(x)U>, (21) 

where a(x) is a smooth complex function of coordinates. Using equations (|Tl]), it's readily seen that 
a cannot be an arbitrary function of coordinates, it rather satisfies the equation 

i A 'V A A>oc = 0, (22) 
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from which in view of 2-spinors' properties follows 

^AA'Ot = pa£a> (23) 

for some p A and, consequently, - the eikonal equation for a{x) 

V AA 'aV AA ,a = 0. (24) 

Before we carry on, we need to establish an auxiliary result. Let us rewrite equation (|i"7| ) in partial 
derivatives 

Vaa>t b = X BB '<S> B , A U> + U'Sa- (25) 

Using the orthogonality £ A '£ A ' = 0, we immediately verify the validity of the following relations: 

V AA 'a B 'V AA ,T B = 0, (26) 
V aa 't b V aa ,t c = 0, (27) 

which along with equation (|9|) lead to the eikonal equation for any function \{T a ) of twistor components 

V AA '\(T a )V AA >\(T a ) = 0. (28) 

Going now back to the main goal of this section and taking into consideration the eikonal equation 
for a(x), we are guided by the equation (|28l ) to conjecture the following 

Preposition. 2 Transformation of the type are symmetries of (j^) iff a is a function of T a and 
&B'A transforms according to 

$B>A -> & B 'A = ®B>A + V AB > In a. (29) 

Replacing in (@) £ A > and <&b'A by their transformed values, after some simplification we obtain the 
following condition of form-invariance of d8|) : 

^B'^AA'Ct - iA'^AB'Ot = 0, (30) 

which is skew-symmetric in A',B' and therefore equivalent to the equation (|22|) for a(x). Taking 



now into account equations (g) and (25) and carrying out simple manipulations we show that if 
a = a(T a ) = a(£ A ',T ) then equation ( [22] ) is identically satisfied. This proves the sufficient part of 
the preposition. 

To prove the converse, suppose that transformation (^Tj) is a symmetry of (||). This yields the 
following: 

i B 'da = a(*BM - ®b>a)w a , (31) 

where, as before, w A = dX AA ' £ A /. Making use of ( |i~6| ) and (0), it's easy to see that the exterior 
product of equation (|3l]) with any two differentials of the twistor components vanishes, leading to the 
functional dependence of a on the corresponding twistor components. More symmetrically, this result 
can be expressed as asserted in the preposition, i.e. a = a(T a ). Passing then to partial derivatives in 
(|3l| ) we obtain the relation 

Cb'Vaa' In a = {&b'A ~ ®b>a)U> , (32) 



7 



from which the desired transformation rule ( p9| ) for the "potentials" &b'A follows (in a special case 
when the primed subscripts are equal). This completes the proof of preposition. 2. ■ 
Some words are in order about the nature of transformations 

= acne- (33) 

which may be called restricted, or weak gauge transformations 28]. We remark that according to its 
definition ( |l4|) the conjugate spinor t b transforms in a similar way 

t^t' = a{T a )r (34) 

so that both (33), (pi) together imply that the gauge symmetries of GSE may be considered as the 



transformations of twistors of the form 

T -> T' = a(T a )T. (35) 



It's clear that composition of transformations ( pq ) is a transformation of the same type. Their 
associativity is also obvious. But the existence of inverse transformation is not so evident. However, 
according to preposition. 1, T a and its image T' a both have only two functionally independent compo- 
nents, and a(T a ) depends, essentially, on these two components (of T a ). So we can always express the 
two independent components of T a through those of T' a . Substituting them in a _1 (T a ) results in the 
inverse transformation T = a -1 (T ,a )T' which is of the same type as (|35|). Hence these transformations 
constitute a groupl In fact it is a proper subgroup of the full C(l)-gauge group of transformations (^Tl), 
the latter itself being generally not a symmetry of the GSE . The statement that this subgroup is a 
proper one, becomes quite evident if we recall that a(T a ) should be subject to the eikonal equation 

©■ 

Finally, we note that under the transformations (|3^) the trace part of the matrix 1-form A = 
Tr(<frdz) = a* AdX AA with components = ^A'A transforms gradient- wise 

A + dlna, (36) 

as the electromagnetic potential 1-form does under the gauge transformations (this may be also seen 
from the expression flj~2"D). In view of the 4- vector properties of $ under Lorentz transformations (|6|) 
we are brought to adopt the interpretation of the 1-form A as the potentials and of the respective 
gauge-invariant 2-form F = dA as the electromagnetic field strengths (of course, up to an arbitrary 
scale factor only). In the following section we obtain Maxwell equations for this 2-form, therefore 
elaborating its electromagnetic interpretation. 

Let us now look at the gauge transformations (|3^) from the viewpoint of the geometry of twistor 
space. The Abelian nature of the transformations studied results in the fact that the ratio of any two 
twistor components T a is invariant under them. Thus, such transformations are projective in origin. 
Not only the planes (|l~4] ) but also the surfaces (18) are form-invariant under transformations (|35| ) and, 



consequently, give rise to another solution of the GSE (with the same electromagnetic 2-form F). So 
we may consider the equivalevce classes of the solutions (and of the surfaces ( |l8|) respectively) which 
can be obtained one from another via the gauge transformations (|35[). That's why we may restrict 
ourselves to consider only the projective twistor space CP 3 . However, projective structure of this 
type differs essentially from the conventional one which originates from the transformations of the full 
gauge group (^l|). We shall return to this problem in section 6, and for the time being shall deal with 
the full structure of the space of (null) twistors. 
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5 Integrability conditions, self-duality and the source-free gauge 
equations 

As was mentioned in the introduction, the GSE (|lj) may be viewed at as the condition that must be 
met for a spinor to be covariantly constant with respect to the IB-connection 1-form 

= <5>dz. (37) 

It may be interesting to note that in the 4- vector representation B-connection (37) turns into the affine 
connection of the form M, ^] 

^ p = 6^„ + 5$* v - ^ - < pA $ A (38) 

which gives rise to the effective complexified geometry of Weyl-Cartan type. For this B-induced 
geometry the Weyl vector of nonmetricity and the pseudotrace of the torsion tensor appear to be 
proportional to each other and are both expressed via the components of the basic gauge field $(x) []. 
According to definition (|37|), the initial GSE (||) may be rewritten as follows: 

d£ = Q£ (39) 

The GSE is overdetermined (8 equations for 6 unknown functions) and both the spinor and the 
"connection" gauge fields are to be determined from it. Dynamics of the connection field Q(x) can be 
obtained by external differentiation of (|39|) which yields 



R£ = (d£l - n A fi)£ = 0, (40) 



where in parentheses the matrix curvature 2- form R of the connection ( |37| ) arises. Since the spinor 
£ is not arbitrary but subject to ( |39"| ) the integrability conditions ( f40| ) don't imply the triviality of 
curvature [], instead they lead to its self-duality Jl|, |2|. 

To demonstrate this, we note that for connection (|37| ) the curvature R is of the following, rather 
specific form 

R= (d$ - <£>dz$) Adz = ir Adz, (41) 
where a new B-valued 1-form tt emerges, with the components 

TTA'C = TT A >CBB'dX BB ' = (V B B'*A'G " <5>AiB$Bic)dX BB ' . (42) 
The integrability conditions (|40| ) take then the form (ir A dz)^ = or, in matrix notation 

n A , C BB>dX BB ' A dX cc 'ia = 0. (43) 
Making use of symmetry properties we obtain from the last relation 

*A' C C(B'£c') = 5 ( 44 ) 



6 Vector fields covariantly constant with respect to the torsion- free Weyl connection have been studied in [^6|; they 
are closely related to the symmetries of Weyl manifolds jl7j|. Relations between the nonmetrical and torsion parts of 
such connections were considered in It. 

7 At this point our approach differs essentially from that of Buchdahl |l8| , Penrose and Plebanski J2^] who assumed 
that the integrability conditions like ( |40[ ) should be satisfied identically for an arbitrary spinor field, in order to ensure 
the existence of an "exact set" of solutions to field equations. 
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so that for any nontrivial solution £(x) it follows 



•K/FgW = VcB'&A? + $B>C$A? = °- (45) 

Decomposing now in a usual way the curvature (^Tj) into the self- and antiself-dual parts it is easy 
to verify that equations (f45| ) are just the conditions for its self-dual part to vanish, whereas the other 
(antiself-dual) one R takes the form 

Rj(BC) = ^%^A'C) ~ ^%^A>C) (46) 

and satisfy additional integrability conditions = (we'll not make use of them below). 

Thus, though the curvature 2-form (|4l| ) of the connection 1-form (|37| ) is not identically (antiself- 
dual (i.e. self-dual in the "strong" sense), it nesessarily becomes (anti) self-dual on the solutions of the 



GSE. For this reason, the property has been called in [27] weak self-duality. 

Physically, the expression (Jt6|) represents the strength of a matrix gauge field; in particular, its trace 
part Fbc = (bc) = ^ A (B^A'C) corresponds to the aforedefined electromagnetic field strength 
F = dA whereas the trace-free part of ( |46| ) defines the strength of a Yang-Mills type ' field 

Indeed, in view of Bianchi identities 

dR = Q A R — R A Q, (47) 

self-duality of curvature R+iR* = leads to the source- free Maxwell equations for the electromagnetic 
2-form F = Tr(R) = R A A' 

dF* = = dF = 0, (48) 

and to the equations of Yang-Mills type for the trace-free part of curvature F A P = R A P — \F5 ® . 

Generally speaking, electromagnetic 2-form F is a C-valued field, yet in view of its self-duality it 
reduces to an R- valued 2-form F related to F by 

F = F — iF*, (49) 

for which Maxwell equations do hold as well, so that the number of its degrees of freedom is just equal 
to that of a usual real electromagnetic field. Explicitly, from the symmetric part of the integrability 
conditions ([45|) we get for the C-valued "electric" E and "magnetic" H field strengths 

E + iH = 0, (50) 

so that %(H) = M(E), = and the pair represents an M-valued electro- 

magnetic field subject to Maxwell equations. Note that from the skew-symmetric part of (^5|) in 
the 4- vector form we get the following " inhomogeneous Lorentz condition" for the C-valued 

electromagnetic potentials = &a'A : 

dpA* + 2A^ = 0, (51) 

which is also identically satisfied on each solution to the GSE. Condition (^) is by no means gauge 
invariant in a usual sense but it is invariant with respect to the weak gauge transformations ( |33| ) 
provided the potentials satisfy the GSE. 

8 Owing to the restricted (weak) gauge symmetry this is not precisely what is usually understood as the YM fields; 
however, the dynamical structure of the gauge equations is completely the same, the restrictions are imposed only on 
the solutions. 
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As to the fields of Yang-Mills type, they may be always expressed via the electromagnetic field 
strengths and the spinor £4/ itself and for this reason can't be regarded independently. Note also that 
the real or imaginary part of the trace-free curvature F A ^ being taken separately would not satisfy 
the source-free YM equations, in view of nonlinearity of the latters. Therefore, the YM-like fields 
here are necessarily complex. Other details about the pecularities of YM fields in the framework of 
algebro dynamical approach can be found in 



6 GSE and the null shear-free geodesic congruences 

Let us recall that via the elimination of potentials &a'A the GSE (|l]) takes the form (|Tl|). Once its 
solution £a'( x ) is found, a field of a null 4-vector ku(x), k^k^ = can be defined as 

k = k^dx^ = UU'dX AA ' . (52) 



Vector lines of this field define a null congruence for which the shear- free criterion [14] 

Z A 'Z B 'Vaa<b> = (53) 
follows readily from (|Tl|). Hence each solution £a' °f the GSE actually defines a null shear-free 



geodesic congruence (SFC). Contrary to (11), the SFC equations ( p3[ ) are invariant under the full 
complex Abelian gauge group (^) and reduce to the system of two equations in partial derivatives 

V W G = GV U G, V V G = GV^G, (54) 

where G = £i'/£o' is the gauge invariant, and the following generally accepted notation has been used: 

r AA' ( u w \ _ f x° + x s x 1 — ix 2 



X AA ={_ " = * 2 3 (55) 
\w v J yx + ix x — x I 

Four real quantities {x*;x }, i = 1,2,3 correspond to the Cartesian space and time coordinates 
respectively. Note that as for the individual spinor components ^a'i they remain fully indeterminate 
by the SFC equations which impose restrictions only on their quotient G{x). 

Let us compare this with the GSE (|Tl|). The latter is equivalent to the system of four equations 
for only two spinor compoments £4/ 

VwU> = GVuW, VvU> = GVuU (56) 



from which again the equations (M) follow for the quotient G(x). Multiple solutions of ( |5q ) with the 
same G correspond to different potentials but have the same strengths of the electromagnetic and YM 
fields assosiated. In view of this, further on we identify the GSE and SFC equations by considering 
only the projectively invariant part of the GSE represented by system ( |54[ ) (one may regard this as a 
fixing of the gauge £0' = !)• 

General analytical solution of (|54] ) for G(x) immediately follows now from the preposition. 1. in 
the form of an algebraic equation 

H(G, t°, r 1 ) = H(G, u + wG,w + vG) = 0, (57) 

which determines implicitly the function G(x). Here n is an arbitrary holomorphic function of three 
complex variables. Equation ([5?]) manifests the functional dependence of three components G, r°, r 1 of 
a projective twistor T a assosiated with the solutions of GSE. For the SFC equations (53) the equivalent 
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result is well known as the Kerr theorem [14]. Note that the solutions of (^4|) in the form ( p7\) are 
defined exept at the points of the singular set whose equation ((2^) simplifies now to 



By multiplying the two equations of (p4J) we obtain once more the 4-eikonal equation for G(x) in 
the form 

V U GV V G - V w GVu,G = 0, (59) 
while by differentiating them we verify that G(x) satisfies also the linear d'Alembert equation 

UG(x) = (V U V„ - V w Vu,)G(x) = 0. (60) 
Note that in view of (|59|) every C 2 -function A(G) is also harmonic on the solutions of the GSE, 

□A(G) = 0. (61) 



Using now ansatz ( |12| ) for the potentials <&a'A an d taking into account (|5^), we can express the 
electromagnetic field strengths via the 2-nd order derivatives of In G as 



T)0 



V^V^lnG, F X i = V v V w lnG, F 01 = V W V^ InG, (62) 



so that satisfaction of the source-free Maxwell equations is then ensured for ( ]62| ) in consequence of the 
wave equation (|6l| ) for A = InG. Differentiating twice the identity (]57| ) with respect to the (spinor) 
space-time coordinates we obtain finally for the strengths (|62|) the following symmetric expression: 



where the function P is defined by fl58|) and {n^, n^e}, A, B = 0, 1 denote the (1-st and 2-nd order) 
partial derivatives of II with respect to its twistor arguments r°,r . We return to this expression 
below. 

Close relation of the GSE to the SFC makes it possible to introduce one more geometrophysical 
structure - an effective Riemann metric. In fact, it's well-known [31, |3i| that we can transform the 
flat Minkowski metric rj^ u into a metric g^ u of the Kerr-Schild type 

= Vnv + hk^k u , (64) 

and that the main characteristics of SFC (geodesity, twist and shear) are preserved under this 
transformation. Here h is a scalar field and the congruence k given by ( |52"|) acquires a projectively 
invariant form 

k = du + Gdw + Gdw + GGdv, (65) 



G being complex conjugated of G. Now we resort to the results of the classical paper plj] where it has 
been proved that the metric (|64|) satisfies Einstein-Maxwell electrovacuum system for any function G 
obeying algebraic constraint (|57|), with the function II linear in twistor arguments r ^ 1 : 



n = <p + (qG + s)r l - (pG + q)r°. (66) 

Here 92 = (p(G) is an arbitrary analytic function of the complex variable G, s and p are real constants 
and q is a complex constant. Not going into details we just note that according to the results of [31] 
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the scalar field h in ( |64| ) is determined, up to an arbitrary real constant, by the function IT and some 
another function ^(G) independent of <f(G) and related to the electromagnetic field arising therein. 
These electromagnetic fields are defined in the curved space-time with metric ( 34 ) , and generally they 
are different from those emerging in our approach and satisfying Maxwell equations in the flat space- 
time]^]. However, for the most fundamental Reisner-Nordstrem and Kerr-Newman solutions these fields 
coincide, the only difference being in that in our approach the electric charge is fixed in magnitude by 
the GSE itself (see the next section). 
In g|, it; 

was shown that singularities of Riemann curvature of the Kerr-Schild metric (p4|) 



correspond just to the condition (|58|). On the other hand, the expression (|63|) demonstrates that the 
same condition P = determines the set of singular points of the electromagnetic field. It may be 
verified that this is true also for the strengths of Yang-Mills fields assosiated with the solutions of the 
GSE. 

Hence, with each solution of the GSE an electromagnetic, a C-valued Yang-Mills and an effective 
gravitational field can be naturally associated, their singularities are all determined by equation ( |58|) 
and therefore coincide in space and time. This makes it possible, in the framework of the algebro- 
dynamical approach based on the GSE, to consider particles as common singularities of all fields 
involved. This general concept will be developed in the next section. 



7 Quantization of electric charge and "particle-like" solutions of 
GSE 

We'll briefly review now the main solutions of the GSE known to date which all may be obtained by an 
appropriate choice of the function n and solving subsequently the algebraic equation ( |54] ) . In order to 
find the solutions in a simple explicit form one usually restricts the consideration to functions n which 
are quadratic in G (linear functions result in zero electromagnetic strengths (^)). The fundamental 
static solution is generated by a function 

n = Gt° - t 1 + 2ia = G(u + wG) - (w + vG) + 2ia = 0, 

(a = Const G M), from where we get 

G = - ( % ee fi±M (67) 

(z + ia) ± r* ( z + ia) ± x 2 + y 2 + (z + ia) 2 



Electromagnetic field (32) assosiated with this solution 



E-iH = ±-^; (E + iH = 0), 

where = {x, y,z + ia} has a ring singularity of radius a, an only possible electric charge q = ±1/4 
(in the dimensionless units we use) , a dipole magnetic and a quadrupole electric moments equals to qa 
and qa 2 respectively [^] . Apart from the restriction on charge, the electromagnetic field ( |68| ) together 
with the Riemann metric assosiated with (^) via the SFC (pBD, accurately reproduce the field and 
metric of the Kerr-Newman solution (in the coordinates used in |3l|). Particularly, for a = the 
solution ( |67[ ) represents the stereographic map S 2 — » C whereas the fields turn into the Coulomb one 
and the Reisner-Nordstrem metric respectively. 

9 At the same time they both are generally different from the fields which may be defined for the SFC using the 
Penrose twistor transform, see |l4], EH| . 
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Quantization of electric charge seems to be a profound property of the GSE-solutions discovered 
in ||]. It is a consequence of self-duality condition (p5(i| ) which together with the gauge symmetry 
of GSE ensures the relation q = N/4, N € Z for the values of electric charge assosiated with every 
solution of the GSE. This property has both topological and dynamical origins, the latter being related 
to the overdetermined structure of the GSE. The proof of the general theorem on charge quantization 
will be presented elsewhere. Contrary to the recently developed |34|, |35| approaches to the problem 
of quantization of electric charge, which are purely topological, in the framework of the GSE the 
charge of the fundamental static solution (|37j) can be of only one fixed and minimal possible value 
and, therefore, can be naturally identified with the elementary charge. Together with the well-known 
property of Kerr-Newman solution to fix the gyromagnetic ratio g = 2 equal to the ratio of the 



Dirac particle [36] the appearance of elementary electric charge within the theory makes it much more 



legitimate to interprete the fundamental solution (67) as a classical model of electron ( in comparison, 
say, with the models of Lopez 37], Israel [38] and Burinskii [39] based on Einstein-Maxwell theory 
itself). 

According to a general theorem proved in |52] static solutions of the SFC (and thus of the GSE) for 
which the singular set is bounded in 3-space (below we call them particle-like [pQ[ ) are all exhausted by 
the Kerr solution (|67|) (up to 3-translations and 3-rotations) . If, however, we remove the condition for 



a solution to be static and get out of the class of functions (pq) dealt with in [31] we discover a lot of 
time-dependent " particle- like" solutions with bounded singularities of different dimensions, 3-shapes 
and time evolution. 

In particular, an axisymmetric solution of a pacticle-like type generated by the function 

n = tV + b 2 G 2 = 0, b = Const, (69) 

has been found in [^ , p9f [. For real b it corresponds to the case of two singularities with elementary 
charges +1/4 and —1/4 undergoing head-on hyperbolic motion for which the electromagnetic field 

^- ± ^s7^' ^- ± ^s7^' ( 7 °) 

is identical to that of the Born solution. Here the following notation is used: 

p 2 = x 2 + y 2 , s 2 = t 2 -z 2 , M = s 2 + p 2 + b 2 , A = M 2 — 4s 2 p 2 , 

and singularities are defined by the condition A = 0. For imaginary b = ia, a£l one has at t = a 
neutral ring- like singularity which then expands to a torus. After an interval of time t > \a\ singularity 
turns into a self-intersecting torus depicted in Fig.l. 



We mention also a particle-like solution for which the singularity has a plane 8-figure shape at 
t = 0, as well as a wave-type solution with a helix-like singularity |3(J. The latter stands for the 
analogue of electromagnetic waves in electrodynamics based on the GSE. 

A beautiful representation of the solutions of SFC-equations has been proposed by E.T.Newman f4(i[| 
and developed then in the works A. Burinskii, R.P.Kerr and Z.Perjes |41[| . They regarded a point-like 
source moving along some curve in a complexified Minkowski space-time CM. Then a trace of its 
complex light cone on the real space-time M is just a null congruence which turns out to be shear- free. 
Kerr congruence is the simplest example of this representation (a point-like source resting at some 
place of the orthogonal to M imaginary subspace of CM). The examples above-presented demon- 
strate, however, that the structure of singular set for such "complexified Lienard-Wiechert field" can 
be quite nontrivial. 
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t=0 t> I a I 




Figure 1: Singular set of electromagnetic field (ff(]) for neutral solution ( |69[ ) at initial (t = 0) and final 

(t > \a\) stages 

To conclude, a number of exact solutions of the source-free Maxwell equations may be obtained 
in a purely algebraic way, some of them being unknown before. They are defined exept at the points 
where the electromagnetic field turn to infinity. These points constitute a set which may be 0-, 1- 
or 2-dimensional and for a particle-like solution is bounded in 3-space. Generally, it's impossible 
to cover such set by some 5-like source (because of the multivaluedness of the Kerr- type solutions). 
Nonetheless, the quantum numbers, 3-shape and time evolution of the singularities are well defined 
and nontrivial owing to the "hidden nonlinearity" of Maxwell equations in this theory, i.e. to their 
origination from the primary nonlinear GSE. The latter ensures some "selection rules" to exist for the 
solutions of Maxwell equation, in particular the restrictions for the allowed values of electric charge 
and the violation of the superposition principle (generally, the superposed solution may satisfy linear 
Maxwell equations but not the GSE itself). A detailed discussion of the status of singular particle-like 
solutions may be found in |3C[| . 



8 Conclusion 

In this paper we didn't claim to present an alternative field model or an algebraic method to derive the 
solutions of classical field theory' conventional equations. We only attempted to study the properties 
of differentiable functions of IB-variable themselves, i.e. to construct a new version of noncommutative 



analysis. In the particular case considered, the generic conditions of B-differentiability |l], ||, 24] 
reduce to the GSE (|l|) which naturally involves the gauge and the 2-spinor structures and manifests 
wonderful correlations with the structures and the language generally accepted in classical field theory. 

In fact, we make only three fundamental assumptions in order to physically interpret the abstract 
mathematical scheme we develop: 

1) about space-time as a (4-dimensional subset of the) vector space of IB-algebra, 

2) about physical fields as differentiable functions of IB-variable and 

3) about particles as bounded singularities of strengths (curvatures) of gauge and metric fields naturally 
assosiated with original B-differentiable functions-fields. 

From the physical viewpoint, the GSE may be regarded as some peculiar (nonlinear, non-Lagrangian, 
overdetermined) field model dealing with effectively interacting 2-spinor and electromagnetic fields, 
the dynamical equations for the latter being not postulated but derived from the GSE itself as its 
integr ability conditions. 

Twistor structure arises in the theory in a natural "dynamical" way via the integration of the GSE 
and makes it possible to obtain all its solutions, as well as correspondent solutions of gauge equations, 
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in a fairly simple algebraic way p°| 

Explicitly, by the examination of algebraic equation (]57| ) a wide class of exact solutions, of linear 
Maxwell equations in particular, may be found, solutions with bounded singularities among them. 
Condition (^) performs the role of equation of motion for this particle-like objects but at the same 
time predetermines the charactetristics and the spatial shape of singularities realizing in such a way 
Einstein's hypothesis of super- causality J23|| . 

Due to violation of the superposition principle for the original GSE, the evolution of the particle- 
like objects simulates physical interaction, and the dynamical perestroikas of the structure of singular 
sets may be interpreted as the transmutations of particles. All this reveals close relations to the theory 
of singularities of differentiable mappings and to the catastrophe theory plfl . For example, in GTR the 
singularity condition (^) is recognized as the condition of caustics of corresponding light-like beams 
formed by the shear- free congruences (|65|). 

It seems also that at least some remarkable properties of the GSE may have consequences for field 
theory in general. In particular, we refer here to 

1) a possible expansion of a class of gauge models in account of weak gauge symmetry (^) discovered 
and using the connections of the form (^) , (^8|) ; 

2) a natural possibility to obtain the selection rules for electric charge, spin etc. using the overdeter- 
mined structure of field equations; 

3) a total algebraization of the primary PDE equipped by a twistor structure using the analogue of 
Kerr theorem or its generalizations like (|i~8|); 

4) a possibility to establish the structure and evolution of singularities without even explicitly ob- 
taining the solution of field equations itself (via elimination of the main field function G(x) from the 
system of algebraic equations ( |57| ) and (p0|), a procedure proposed in p2f). 

There are at least three ways to consider the material presented in the paper and the role of 
GSE in particular: as a beatiful mathematical toy, as a powerful method to generate the solutions of 
conventional field equations, and as a fundamental dynamical system (or an example of such) primary 
to conventional Lagrangian systems. The utilization of classical dynamics based on the overdetermined 
structures like GSE requires quite new methods of quantization. Alternatively, one can claim to explain 
the quantum properties on the whole via, say, the stochastic behaviour of an ensemble of particle-like 
(singular) solutions or by other means of a consistently classical consideration. 

In any case, to find the correct approach to quantization and to physical interpretation in general 
one needs to study accurately the properties of the classical solutions themselves [^]: their complete 
classification, dynamics and bifurcations. All these problems are evidently related to the general 
theory of singularities of differentiable maps [21, p2j |. Nevertheless, the already discovered properties 
of B-differentiable functions- fields and of numerous geometrophysical structures they give rise to, look 
like somewhat striking and bring one back to Pithagorean philosophy about the numerical origin of 
fundamental physical laws. 



9 Appendix 

The below reviewed approach to differentiability in the quaternion-like algebras have been motivated 
by the old works of G.Sheffers |ll| (see also |l^]) on the analysis over an arbitrary commutative 
associative algebra A, and is a direct generalization of Sheffer's approach to the noncommutative case. 
Let F(Z) be an A-valued function F : A i— ► A of an A-variable Z € A. Sheffers defines the condition 

10 In R.Penrose's twistor approach we need to integrate a function of twistor variable in order to get a solution to the 
wave equations; contrary, in our approach even this is not necessary. 

11 To clarify the correspondence with quantum theory the particle-like multisingular solutions are especially interesting. 
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of its differentiability in A using proportionality of the linear parts of increments (differentials) dZ, dF 
of independent variable and its function respectively as 



dF = H(Z) * dZ, (71) 

where H E A and (*) denotes multuplication in B. For division algebras ( |7l| ) is equivalent to the 
condition of existence and path-independence of the derivative H(Z) = dF * dZ^ 1 = F'(Z), and 
in a particular case of complex algebra C leads to Cauchy-Riemann equations. Generally, however, 
(JTT|) may be applied to the algebras with zero divisors, in particular to that of the double and dual 



numbers. Linear differential equations relating the components of F{Z) follow from (71) via the 
elimination of H(Z) and fully resemble the CR-equations for the functions of complex variable. In 
many aspects the analysis constructed by Sheffers is quite similar to the complex one, so that a wide 
class of A-differentiable functions subject to (|7l| ) can be found, including all the polinoms in particular. 

In the r.h.s of (fn]) an invariant A-valued differential 1-form of the most general type is present, 
which can be constructed via the algebraic operations in A only. A natural generalization of (|7l|) to 
the case when A is noncommutative (yet associative) seems to be the following condition (see U 
and the references therein): 

dF = L(Z) *dZ* R{Z) (72) 

for the function F{Z) to be differentiable in A. Here L,R € A are the so called left and right semi- 
derivatives of F(Z) respectively. For a given function F they are determined (if exist) not uniquely, 
but at least up to a transformation L — > aL, R — ► a~ 1 R, where the function a(Z) takes the values in 
the centre (the commutative subalgebra) of A. 

For commutative algebras ([72]) reduces back to ( |7l| ) with H(Z) = L(Z)*R(Z) . On the other hand, 
if in noncommutative case we take, say, R = e (we assume the unit element e to exist in A), we come 
to the condition ( |7l| ) with H(Z) = L(Z). Nonetheless, at least for the algebras of quaternion type 
this condition is known to be too restrictive, being satisfied only by the linear function F = A* Z + B, 
where A, B are constant elements of algebra (see e.g. || |l0|1). 

As to the general A-differentiability condition (|7^), it defines a wider class of functions. Partic- 
ularly, for Hamilton quaternions H condition ([72] ) turns out to be an algebraic analogue of for the 
mapping Z i— > F(Z) to be conformal in E 4 p5| p6| , In this regard equation ( |72| ) can be viewed 
at as a natural generalization of complex holomorphy. However, in E 4 conformal mappings constitute 
only a finite 15-parameter group, contrary to the infinite-dimensional complex case. Thus, for division 
algebra EI the class of H-differentiable functions, as defined by (|72|), is again too narrow to be used, 
say, in field dynamics. 

The situation changes radically when we come to consider noncommutative algebras with zero 
divisors, in particular the algebra of biquaternions B (quaternions over C). For simplicity let us below 
restrict ourselves by consideration of the full N x N matrix algebras over M or C (for N = 2 we have 
just an isomorphism Mat(2, C) = B). Then for a determinant of the matrix of differentials dF in the 
l.h.s. of ( |72"P we obtain 

det \\dF\\ = det \\L(Z) * R{Z)\\ det \\dZ\\ = X(Z) det \\dZ\\. (73) 



In the case both L,R are invertible, so that X(Z) ^ 0, relation (|73| ) defines a conformal mapping 
with a scale factor \(Z) and a (positively indefinite or complex) infinitesimal metric interval repre- 
sented by corresponding determinants in (|73|). Particularly, for B we deal with conformal mappings 
in complexified Minkowski space CM. 

In a remarkable way, however, for det L = or similarly det R = 0, we have X(Z) = and relation 
( |73| ) defines a reduction of the full vector space of A to the subspace of null elements (to the complex 
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"light cone" for B). Such mappings may be called degenerate conformal mappings; they constitute a 
wide and important class: in the context of the presented theory just these mappings ( differ entiable 
A-functions) are identified with the physical fields. 

In the N x N matrix notation (fF3) takes the form (A,B, ... = 1, ...N) 

VabFcd = LcaRbd (74) 

where V ab stands for a derivative operator with respect to the coordinate Z AB . For some indices 
C, D being fixed we denote Fqd = ^, Lqa = <f>Ai Rbd = ipB, and the relation ([74]) becomes 



V A b £ = <t>AipB (75) 

In view of the zero determinant of the matrix in the r.h.s., we get the equation 

det ||Vab£|| = 0, (76) 

which have to be satisfied for each matrix component Fcd = S E R, C of an A-differentiable function. 
Equation (76) is a nonlinear analogue of the Laplace equation in the complex analysis, and nonlin- 
earity arises here as a direct consequence of the account of noncommutativity in the definition of a 
differentiable function (|72[). For the case of biquaternions B (76) represents a (complexified) A-eikonal 
equation 

(V 00 S)(VnE) - (V iS)(V 10 S) = (77) 
which in Cartesian complex coordinates z°, z 3 = z 00 ± z , z 1 , z 2 = z 01 ± iz 10 takes a familiar form 

<!§> 2 -<§> 2 -<II> 2 -<f§> 2 ^ < 78 > 

In the paper we restrict our consideration by a particular, yet a basic subclass of A-differentiable 
functions for which one of semi-derivatives, say R(Z), is proportional to the function F(Z) itself. 
Redefining then L(Z) — ► $(Z) we get instead of (f72|) 



dF = $(Z)*dZ*F(Z). (79) 

Let now {^ c ^}, C = 1, ...N be N columns of the matrix F(Z); then we can present (7£) in a form 
of a system of matrix equations 

= QdZtW (80) 

(here and below we omit the symbol of matrix multiplication), which all have to be satisfied with the 
same (left semi-derivative) matrix <&(Z). The quantities £(Z) are evidently SL(N, C)-spinors with 
respect to the symmetry transformations of ( |30| ) (for JV = 2we have demonstrated this in section 2). 

Different spinors } may be either functionally dependent or not 0: in any case an arbitrary 
solution of ( |79| ) may be constructed from (and decomposed into) a set of N solutions {<J>(Z), ) {%)} 
of the system 

d£ = §dz£. (81) 

Conversely, from a solution to ( |8ll) we easily obtain at least one class of the solutions to the original 
system ( |79| ) by setting N spinors ^ C \Z) to be globally proportional to each other (or zero, exept for 

12 In the latter case we return back to the conformal maps in the JV 2 vector space of A, with the norm represented by 
the determinant. 
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one of them) . Eliminating the semi-derivative matrix from the overdetermined system ( p\\ ) we come 
to a nonlinear system for the components of £ (i.e. of B-differentiable function) only which resembles 
in a certain sense the Cauchy-Riemann conditions in the complex calculus (for N = 2 this is just the 
system ( |Tl| ) of the paper). 

Thus, we have shown that (in a particular case R{Z) = F{Z) ) the A- differentiable functions, 
as defined by ( |72| ) , can all be found from the condition flBl] ) . For the algebra of biquaternions B this 
equation becomes equivalent to the GSE (H) studied in the paper if we only assume the coordinates 
to be real-valued, so that the matrix Z —* X, X = X + is considered to be hermitian. 

This is the only ad hoc assumption which is motivated by physical considerations and which does 
not follow from the algebraic structure or the differentiability conditions ( f72|) themselves (precisely, 
we have to deal with the full structure of 4-dimensional complex space). Under the assumption made, 
the coordinate space reduces to the Minkowski one and the whole theory, including the fundamental 
4-eikonal equation (|78|), becomes Lorentz invariant. 

Some details and generalizations of the approach afore-presented can be found in the monograph 



|IJ and, partly in English, in |25|, |, |4|. 
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